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ABSTRACT 

This  thesis  presents  a  discussion  of  the  problems 
involved  in  estimation  of  Circular  Error  Probable  (CEP) , 
Several  estimators  are  compared  through  simulation,  under 
two  models;  the  power  model  and  the  Rayleigh  model.   Several 
measures  of  effectiveness  are  calculated  for  each  of  the 
competing  estimators.   It  is  found  that  maximum  likelihood 
estimation  based  on  the  power  distribution  performs  well 
for  "heavy  tailed"  distributions;  the  Rayleigh  unbiased 
estimator  performs  well  for  most  other  situations. 
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I.   INTRODUCTION 

Evaluation  of  a  weapons  system  is  frequently  an  expen- 
sive and  time  consuming  task.   Only  relatively  few  tests 
may  be  allowed  because  of  the  high  cost  of  testing.   It  is 
possible  that  a  good  weapon  system  could  be  rejected  because 
of  inefficient  utilization  of  the  small  amount  of  data 
available.   Hence,  efficient  use  of  the  limited  data  is 
required  to  reduce  this  risk. 

Several  authors  have  considered  the  measurement  of 
weapon  delivery  accuracy  using  Circular  Error  Probable 
(CEP) .   This  measure  of  accuracy  is  defined  to  be  the  radius 
of  a  circle  centered  at  the  target  which  on  the  average  con- 
tains fifty  percent  of  impact  points,  that  is,  the  median 
radial  miss  distance.   A  functional  relationship  between 
CEP  and  the  distribution  parameters  is  given  by  the 
following  integral: 


CEP 
F(CEP)  =   /    f(r)dr  =  0.5  , 
0 


where  f(r)  is  the  probability  density  of  radial  miss 
distance . 

Various  parametric  models  for  the  distribution  of  impact 
points  about  the  target  have  been  used  in  connection  with 
such  estimations.   A  discussion  of  these  models  is  given  by 
Jordan  [1] .   But,  deficiencies  in  the  models  have  been 


noted.   Specifically,  they  tend  to  provide  poor  fit  to  the 
"upper  tail"  of  the  radial  distributions.   That  is,  they 
provide  less  than  satisfactory  explanation  of  impacts  ob- 
served far  from  the  target.   A  commonly  used  model  is  the 
"Rayleigh  model",  in  which  impacts  in  the  target  plane  are 
assumed  to  be  distributed  as  the  unbiased  circular  normal 
distribution.   One  might  anticipate  that  Rayleigh-based 
CEP  estimators  might  not  perform  adequately  with  a  large 
proportion  of  actual  test  data.   Several  potential  competitors 
to  Rayleigh  CEP  estimators  were  therefore  developed  for 
comparison. 

Barr  and  Jordan  [2]  developed  and  proposed  the  power 
distribution  m.odel  which  is  based  on  an  infinite  mixture  of 
Rayleigh  radial  miss  distance  distributions.   This  distri- 
bution has  properties  compatible  with  the  interpretation  of 
a  in  the  underlying  Rayleigh  distribution  model;  they  sug- 
gested the  use  of  maximum  likelihood  estimates  of  the  param- 
eters of  the  power  distribution  which  in  turn  provide  a 
maximum  likelihood  estimate  (MLE)  of  CEP. 

The  jackknife  is  a  technique  which  v;as  proposed  by 
Quenouille  [10]  for  reducing  bias  in  estimation  problems. 
We  tried  the  jackknife  procedure  to  estimate  CEP  under  the 
Rayleigh  model. 

The  sample  median  was  also  used  to  estimate  CEP.   Since 
CEP  is  defined  as  the  median  radial  miss  distance,  a  natural 
estimator  for  CEP  is  the  sample  median. 


To  measure  the  performance  of  the  estimators,  a  variety 
of  data  sets  from  known  distributions,  such  as  uniform, 
Rayleigh,  pov;er  and  Weibull  were  generated  by  simulation. 
Samples  from  the  known  distributions  were  generated  to 
test  the  capability  of  the  estimators  to  handle  widely 
disparate  situations. 


II.   MAXIMUM  LIKELIHOOD  ESTIM/iTION  BASED 
ON  THE  POVJER  DISTRIBUTION 

Barr  and  Jordan  [1]/  [2]  have  given  a  discussion  of 
CEP  accuracy  measurement  under  a  pov;er  distribution  model. 
They  showed  that  estimating  the  parameters  of  the  power 
distribution  by  the  method  of  moments  may  not  be  adequate; 

moreover  it  is  possible  only  when  the  ratio  of  the  squared 

—2    2 
sample  mean  and  variance  r  /s    does  not  exceed  7t/4  -  it. 

Consequently,  they  suggested  use  of  the  maximum  likelihood 

(M.L.)  method. 

Analytically  closed  expressions  for  the  M.L.  estim>ators 

cannot  be  obtained,  the  solution  of  normal  equations  yielding 

the  parameter  estimates  requires  numerical  techniques.   In 

this  section  we  examine  a  numerical  method  to  obtain  the 

M.L.  estimates  of  the  parameters  of  the  power  distribution, 

which  in  turn  provide  a  M.L.  estimate  of  CEP.   These  estimation 

procedures  require  use  of  a  computer. 

A.   MATHEMATICAL  FORMULATION 

The  probability  density  function  of  the  radial  miss 
distance  R  having  a  power  distribution  with  parameters  k  and 
C   is 


f  (r;K,c)  =  ^  ( ^)'^'^^  '   r  >  0,  k,c  >  0 
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A  C.D.F.  for  radial  miss  distance  is 


Fj^(r;K,c)  =  1  -  ( ^—5-)'^  '•   r  ^  0,   k  ,  c  >  0   (1) 


C  +  r 

The  CEP  is  given  by 


1       1 

CEP  =  Fj^"^(^)  =  [(2'^  -  1)^2 


Because  radial  miss  distance  is  always  non-negative/ 
one  can  equivalently  consider  squared  radial  miss  distance. 

In  the  case  of  the  power  model,  this  leads  to  some  notational 

2 

simplification.   Let  Y  denote  the  squared  miss  distance  R  . 

Then  from  (1) , 


FY(y)  =  1  -  (1  +  ^)"'^       ;   y  >  0  ,  K,c  >  0 

(2) 

fyCy)  =  ^(1  +  l)""""^       ;  y  1  0  ,  K,c  >  0 


Here  ^  and  k  are  parameters  to  be  estimated  from  observed 
impact  data.   The  relationship  between  CEP  and  these 
parameters  is 


CEP^  =  F^  ^(i)  =  [(2''  -  1)U 


The  method  of  maximum  likelihood  is  based  upon  the 
likelihood  function  (|)(y)  which  is  defined  as  the  joint 
density  function  of  a  sample  of  size  n;  i.e., 


11 


n  n   n      Y  ■  ■, 

4)(y)   =    TT  f   (y  )   =  ^    TT   (1  +  -i)"^"-^ 

i=l     ^     c   i=l      ^ 


Since  it  is  easier  to  deal  with  sums  rather  than  products 
when  maximizing,  <p  iy)    is  transformed  into 

L(y)  =  log  4)  (y) 

n 
=  n  log  K  +  nxlog  ^  -  (k  +  1)  Y.    log(c+y-) 

i=l       ^ 

There  is  no  loss  of  generality  here  because  the  m.aximum  of 
a  positive  function  occurs  at  the  sam.e  point  as  the  maximum 
of  the  logarithm  of  the  function.   The  normal  equations  are 
then 


(K,i;)  =  1^^^  =  ?  +  nlog  c  -  I    log{c+y-)  =  0 


,<  .  .^^ 


and 


dL(y)  _  nK 


g«,0  =  ff^  =  '-f  -   «+i)  .E  ^^-y.  =  0 

^  ^  1=1        -^  1 


The  solution  of  m(K,C)  =  0  and  g{K,c,)    =    0  gives  the 
parameter  M.L.  estimators  k,  C.   As  mentioned  before,  in 
the  present  case,  this  requires  numerical  techniques. 

B.   NUMERICAL  SOLUTIONS:   NEWTON'S  METHOD 

This  method  is  a  Newton's  iteration  technique  and  the 
concept  and  procedure  are  discussed  in  [3].   Newton's  method 
provides  a  statistical  criterion  for  stopping  the  iterations 
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The  pure  form  of  Newton's  method  is  given  by  the 
following  iteration  scheme: 


^j  +  1  "  ^'j  "  tF(X^)]  -"  Vf  (Xj)^ 


where 


X,  =  [K,,L],  f(x,)  =  [iLixL^ikiX).] 


D      D  ^3 


3k     8c 


First  we  compute 


^ii^  =  &+  n  log  c"  -   Z  log  (c^y.) 
3k      k  i=l 


8C       C  1=1  c  +  Yi 


3k  k 

3^L(y)  ^  n  _   ^     1 

3k3c  c    i=l  C  +  y^ 

3c  5        1=1  (c  +  y^) 


Let 


_  _  3  L(y)  .  3^L(y)    ,  3^L(y)  .  2 
3k        3c         3k3c 
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If  D  were  identically  equal  to  zero,  this  method  would  not 
be  applicable.   Now  take 


8^L(y) 

ar 

3^L(y) 
8k8C 

a^L(y) 

9^L(y) 

-1   " 

8k3c 

8k 

D 


In  order  to  develop  the  termination  criterion,  consider 
the  statistic  for  (for  reference  see  [4]) 


B  =  Vf  (X.)F~-^(X.)  Vf  (X.)~ 

3     :    y 


B  is  asymptotically  a  Chi-square  r.  v.  with  2  degrees  of 
freedom,  therefore  it  might  be  reasonable  to  continue  the 
iteration  until  the  statistic  B  becomes  less  than  the  100a 

percentile  point  of  the  chi-square  distribution  with  two 

2 

degrees  of  freedom,  x  (ot)  . 
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III.   MAXIMUM  LIKELIHOOD  ESTIMATION  BASED 
ON  THE  RAYELIGH  DISTRIBUTION  MODEL 

Suppose  impacts,  measured  in  a  suitable  plane,  have 
range  (X)  and  deflection  (Y)  components  of  impact  that  are 

jointly  uncorrected  bivariate  normal  distributed  with  mean 

2 
at  the  target  and  common  unknown  variance  o   in  both  direc- 
tions.  That  is,  for  the  present  development,  it  will  be 
assumed  that  the  joint  density  of  miss  distances  about  the 
target  is  given  by 


r:,         V    ,^      2,-1  -(x^+y^)/2a^ 

f(x,y)  =  (27:a  )   e     ^      '  ,      -»<x,y<oo. 


It  is  easy  to  derive  the  density  of  radial  miss  distance. 

1/2 
Let  R  =  (X+Y)  '       and  transform  to  polar  coordinates  to 

obtain 

,      2ti   r     2" 

Fp(r)  =  P[R£r]  =  — ^—   /    /   R  e  ^o   dR  de 

27TO     0     0 

2 


=l-e       ;   r>0,   o   >0. 


This  distribution  of  the  radial  miss  distance  is  referred 
to  as  either  the  radial  normal  distribution  or  the  Rayleigh 
distribution. 

Then  the  density  function  of  the  radial  miss  distance, 
fj^(r)  is 
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2 

r 

f_(r)  =  -£.  e  ^^    ;   r  >  0,   a^  >  0  . 


To  find  the  relationship  between  CEP  and  a,    we  calculate 


CEP  =  Fj^'-'-ci)  =  y2  In  2  a  :  1.1774  a  , 

which  is  a  well  known  relation.   Hence,  if  we  want  to  estimate 
CEP  with  this  model,  first  we  might  estimate  the  population 
standard  deviation.   The  maximum  likelihood  estimator  for 
a  is  easily  found  to  be 

1 

2n  .  ,  1      ' 
1=1 

where  r, /  . . . ,  r   is  the  sample  of  observed  radial  miss 
distances.   This  estimator  has  a  slight  bias,  however.   An 
adjusted  unbiased  estimator  for  CEP  is 

CEP  =  1.1774  (/n  -L^ )  J,  (3) 

r(j(2n+l)) 

where  the  quantity  in  parenthesis  is  an  unbiasing  factor 
for  a.   It  has  been  shown  by  Chapman  and  Robbins  [6]  that 
CEP  =  1.1848  o  for  sample  size  20  is  a  best  estimator  of  CEP 
under  the  normality  and  other  assumptions  stated  above. 
However  other  estimators,  while  not  as  efficient  as  the  best 
estimator,  may  have  more  robustness  under  variations  in 
the  assum.ed  distribution  of  radial  miss  distance. 
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The  term  "jackknife"  has  been  used  by  Miller  [9]  and 

Quenouille  [10]  to  describe  a  method  of  modifying  a  biased 

estimator  to  bring  it  closer  to  an  unbiased  condition. 

From  (3)  we  see  the  M.L.  CEP  estimator  has  a  bias  term  of 

order  — .   We  are  thus  led  to  try  to  examine  the  jackknife 

/n 
method  to  reduce  this  bias  in  the  estim.ator.   Brillinger  [7] 

has  provided  a  discussion  of  the  asymptotic  properties  of 

estimators  obtained  with  the  jackknife  method  when  it  is 

applied  to  M.L.E's, 

The  jackknife  procedure  is  based  on  dividing  the  data 

into  groups,  obtaining  estimates  from  combinations  of  the 

groups,  and  then  averaging  these  estim.ates.   We  apply  this 

to  our  problem,  as  follows.   The  maximum  likelihood  estimator 

for  a   IS 


-2    1^     2 

2"  j=l   ^ 


Divide  these  data  into  m  groups,  where  r, ,  ...,  r   is  the 
sample  of  n  observed  radial  miss  distances  and  the  size  of 
each  subgroup  is  k.   First  estimate  the  parameter  a   based 
on  all  data 


-  2    1^2 
o     2n    ,  1 

next  estimate  a   based  on  all  data  except  those  in  the  i 
group 
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-21^2      .    ,   _ 
a   =  -^^ —   z.r.    ;   i  =  l,  z ,    ....m, 
-1    2n  .^.  : 


Then  compute  the  pseudo  values 


'^  2     '^  2         ^2 

a.   =mc    -  (m-l)a     ;   i=l,  2,  ...,m. 
1       o  -1  ' 


Finally,  compute  the  average  of  the  pseudo  values, 

_     Z   a. 
8^  =  1=1   ' 


m 


and  estimate  the  CEP  by 

=  y/2  In  2  a 


CEP  = 
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IV.   THE  SAMPLE  MEDIAN 


Since  CEP  is  defined  as  the  median  radial  miss  distance, 
a  natural  estimator  for  CEP  is  the  sample  median  (m) .   By 
definition;  the  median  is  the  50   percentile  of  the  distri- 
bution, i.e.  , 


P(R  £  m)  >_  0.50    and    P(R  >.  m)  >^  .50 

Since  F(CEP)  =  0.50  and  R  is  assumed  to  have  a  continuous 
distribution. 


-11 
CEP  =  m  =  F  -^(i) 


Suppose  that  the  n  radial  miss  distance  observations  of 
the  sample  are  arranged  in  order  of  magnitude,  ^nw  ^(7)' 


•  • 


,  r,  ..   We  shall  consider  only  the  estimates  of  the 


population  median  defined  as  follows: 

(a)  If  the  sample  size  is  odd,  and  n  ==  2p+l,  the  median 
(m)  is  taken  as  the  (p+1)    value. 

(b)  If  the  sample  size  is  even,  and  n  =  2p,  the  median 
(m)  is  taken  as  the  midpoint  between  the  p   and 
(p+1)    values  or 


7^^(P)  "■  ^(p+1)^  • 
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The  distribution  of  the  median  will  tend  to  be  syminetri- 
cal,  and  should  be  very  nearly  normal  in  almost  all  popula- 
tions.  The  distribution  of  the  median  in  the  case  of  small 
samples  from  a  normal  population,  whether  the  sample  size 
is  odd  or  even,  will  tend  to  be  normal  [Hoji  [11]) .   The 
normal  form  is  approached  very  rapidly.   The  value  of  the 
standard  deviations  of  the  median  for  n  even  and  n  odd 
approach  the  same  limit, 


a      ~  1  2    '    o/   n   =    1.25331  o/Jn 


where  o  is  the  standard  deviation  in  the  population  sampled. 
We  use  the  order  statistics  to  examine  the  sample  median 
under  the  Rayleigh  distribution  model  (see  Inselmann  [12]), 
and  find  its  bias  and  variance  so  as  to  enable  one  to  make 
a  comparison  with  other  estimates.   If  n  is  even,  the  dis- 
tribution of  sample  median  is  given  by  the  following: 


^^^(0.5n)^  =  {n!/(n-[0.5n])  !  ([0.5n]-l)  I}f  (r)F(r^Q^^^^)  ^^•^''^  ^ 


r-i  _  p/        ^  ^n-  [0.5n] 
^^        ^^''[0.5n]^J 


For  convenience  we  let  a  =  0.5n,  and 


C  =  nl/(n-[0.5n]) I ( [0.5n]-l) !  =  (2a) 1/a 1 (a-1) ! 
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so  the  sample  median  is  r  .   Then  we  may  write  the  distri- 
bution of  sample  median  under  the  Rayleigh  distribution  model, 


g(ra) 


^2  ^2   ^_^       ^2 

c[-^  exp(— ^-,5-)  ]  [l-exp(— ^)  ]^  [expi--^-^)  ]^ 


-7 

a 


2a 


2o 


2a 


2  2 

r  r    /   ^  (a+1)  V  1  fT       f      ^    \  1^-1 

=  0-2"  [exp(-  i— 2 — )][1  -  exp( j)  ] 

a  2a  2a 


We  will   now   consider   computation   of    the   mean   of   ra, 


E(ra)    =   c   /      ^[exp(-^-i|i^)]  [l-exp(-  ^)]^   ^  dr 
0        a  2a  2a 


=   c 


00      2 

r 


/      ±_[exp(-   ^-%:^^-i)]    ^Z       (-l)^(^"^)  [exp(-^)]V 
0        a  2a  K=0  ^  2a 


^,-'-    ,    T,K-a-l.        f°   r^           ,      r^(a+l+K),     , 
c    I       (-1)     (    ^    )       -^      — T  exp(-   ^^ 5 ^'-)    dr 

K=0  0        a  2o 


and  make   the    following   transformation 


y  =   r^(a+l+K)/2a^    , 


r   =   o  j2y/(a+l+K)       , 


dy  =    [r(a+l+K)/a    ]    dr 
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To  obtain. 


3 
E(ra)  =  ca/T  Z^  (-1)  ^  (^^^)  /  (a+l+K)  ^  f  ./y   e~^   dy 


K=0  0 


3  r-       ^-1  (-l)^(^K^) 
=   caT{^)'/2'      Z   i^-ry   . 

K=0   (a+l+K)  "^ 


We  can  adjust  the  bias  term  and  using  the  same  transformation 
as  above  we  may  now  find  the  second  moment 

2 

o         °°   3       2,  ,,.   a-1  -  ,.K.a-l>  ,    ,  r  .,K  , 

E(r  2)  =  c   /   £-^xp(-^-^±il)   Z   (-^)  (  K  )t^^P(-— 2^]   ^^ 
^        0    0^       20^^      K=0  ^° 

2  a-1  (-l)^(^-^) 
=  2ca    Z   


K=0    (a+l+K)^ 


The  variances  of  this  sample  median  can  now  be  computed 
using  the  well-known  formula 


Var  (r^)  =  E(r^^)  -  [E(r^)]^  . 
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V.   SIMULATION  I-IETHODOLOGY 

The  previous  sections  have  been  concerned  with  a 
development  of  different  estimators  for  the  median  of  the 
distribution  of  radial  miss  distances  obtained  under  a  set 
of  fixed  impact  conditions.   An  analytical  comparison  is 
difficult  to  perform  for  these  estimators,  so  computer 
simulation  was  used.   This  section  suminarizes  the  different 
models  and  estimators  discussed  in  the  previous  sections, 
and  includes  an  analysis  of  the  results  obtained  from  simu- 
lations.  Although  the  sample  problems  do  not  represent  actual 
weapon  test  results,  an  attem.pt  has  been  made  to  generate 
data  from  a  wide  variety  of  "realistic"  distributions. 
Therefore  this  simulation  analysis  should  show  relationships 
between  our  estimators  of  CEP,  and  indicate  their  relative 
robustness,  and  hence  potential  for  use  in  a  variety  of 
applications.   The  alternative  estimators  are  as  follows: 
CEP  1:   maximum  likelihood  unbiased  Rayleigh  estimator 
CEP  2:   maximum  likelihood  Jackknife  Rayleigh  estimator, 
with  subgroup  sample  size  k  =  1 

A. 

CEP  3:   maximum  likelihood  power  estimator  under  Newton's 
method 

CEP  4:   sample  median  estimator 

For  the  purpose  of  our  simulation  comparison  of  these 
four  CEP  estimators,  we  generated  samples  of  radial  miss 
distances  (or,  in  some  cases,  squared  radial  miss  distance) 
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from  several  parent  distributions:   Rayleigh,  uniform,  power 
and  Weibull  distributions.   Throughout,  v;e  assumed  a  sample 
size  of  20  impact  observations.   The  method  of  generation 
was  by  C.D.F.  inversion  of  generated  uniform  (0,1)  values. 
For  each  parent  distribution  (hereafter  called  a  "case") 
we  replicated  estimation  of  CEP's  500  times.   Summaries  of 
the  sample  distributions  of  500  values  of  5  CEP's  were 
developed. 

Eight  cases  were  considered;  in  each,  parameter  values 
were  used  which  gave  a  true  CEP  of  100.   The  cases  are  as 
follows  (U  is  a  uniform  (0,1)  random  variable  throughout): 

Case  I.   (Rayleigh  distribution) 


F_^(y)  =  1  -  e  ^^^^     with   a^  =  lo'^/21n2 

In 


CEP 


=  J2    In  2  a 


2 

Generated  value:   Y  =  -2a   In  U  =  -14426.95  In  U 


Case  II.   (Uniform  distribution) 


Fj^(y)  =  200  y    '     ^'^^^  ^   ^    (0,200) 


CEP  =  100 


Generated  value  Y  =  200  U  . 
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Case  III.   (Power  distribution] 


F  2(y)  =  1  -  (1  +  Z)-<  '•   y  1  0  ,  c,K  >  0 


4 
With  K  =  1  and  ^  =  10 


CEP  =  [^(2*^  -  1)  ]^ 


Generated  value  Y  =  (U  "^  -  1)  ^  =  ^^^  ^h  ~   D 


Case  IV.   (Translated  Rayleigh  distribution) 


K. 


where  2a^  =  14282.6809 


CEP  =  [100  +  2o^  In  2]^ 


2 
Generated  value  Y  =  -  2o   In  U  +  100 


Cases  V-VIII.   (Weibull  distributions) 


Fj^^(y)  =  1  -  e 


w 

V 

-y 


CEP  =  (v  In  2) 


_1_ 

2w 


Generated  value  Y  in  the  table  below 
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Four  sets  of  parameters  (w,v)  were  selected  to  give 
CEP  =  100  and  to  provide  distributions  with  varying  upper 
tail  thickness.   One  case  (with  w  =  1)  gives  an  exponential 
distribution  with  mean  v,  and  this  corresponds  to  the  Rayleigh 
case.   Case  I  above  may  be  considered  also  as  a  Weibull  case. 
The  parameter  pairs  used,  in  order  of  increasing  upper  tail 
size,  are  given  in  the  following  table. 

1 

Case      w      V        Y  =  (-vlnU)^^ 

1 

V  2       10^/ln2   Y  =  (-144260504.1  In  U) ^ 

2_ 

VI  1.5     10^/ln2   Y  =  (-1442695.04  In  U)  ^ 

I  1.0     10^/ln2   Y  =  -4426.95  In  U 

VII  .5      10^/ln2   Y  =  (-144.27  In  U)  ^ 

VIII  .1      10*'^/ln2  Y  =  (-3.623886  In  U)  "''^ 
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VI.   SIMULATION  RESULTS 

As  mentioned  above,  we  tested  four  CEP  estimators  from 
eight  different  cases.   Major  comparison  measures  include 
bias,  variance  and  mean  square  error.   Others  are  included 
in  the  tables  that  follow. 
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CASE  I.   RAYLEIGH  DISTRIBUTION 


Y  =  -  14426.95  In  U 


MOE 

CEP  1* 

CEP  2 

CEP  3 

CEP  4 

MEAN 

99.93 

99.94 

95.96 

107 

STD  DEV 

11.66 

11.68 

12.22 

38.8 

MSE 

136.03 

136.37 

165.77 

1549 

GED  MEAN 

99.25 

99.25 

95.17 

100 

VARIANCE 

136.03 

137.37 

149.45 

1500 

COEFF  VAR 

0.12 

0.12 

0.13 

0.36 

SKEWNESS 

0.15 

0.14 

0.049 

0.35 

KURTOSIS 

0.10 

0.085 

0.307 

0.34 

RANGE 

71.38 

70.87 

83.60 

252 

MINIMUM 

67.25 

67.29 

53.74 

17.7 

.10  QUANTILE  84.82 

84.80 

80.99 

58.4 

.25 

92.12 

92.09 

87.67 

82.5 

.50 

99.63 

99.63 

95.96 

105 

.75 

107.49 

107.34 

104.00 

132 

.90 

114.83 

114.83 

111.32 

156 

MAXIMUM 

138.63 

138.16 

137.34 

269 

CEP  1:   maximum  likelihood  unbiased  Rayleigh  estimator. 

As 

CEP  2:   maximum  likelihood  Jackknife  Rayleigh  estimator, 
with  subgroup  sample  size  k  =  1. 

CEP  3:   maximum  likelihood  power  estimator  under 
Newton's  Second  method. 

CEP  4:   sample  median  estimator. 
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CASE  II.   UNIFORM  DISTRIBUTION 


Y 

=  (200  U)^ 

/\ 

y\ 

/^ 

/\ 

MOE 

CEP  1 

CEP  2 

CEP  3 

CEP  4 

MEAN 

98.61 

96.52 

91.24 

102 

STD  DEV 

9.87 

9.73 

17.08 

41.7 

MSE 

109.32 

107.03 

368.63 

1744.0 

GED  MEAN 

96.09 

96.01 

* 

91.5 

VARIANCE 

97.83 

94.92 

291.89 

1740 

COEFF  VAR 

0.10 

0.10 

0.19 

0.40 

SKEWNESS 

-   0.21 

-  0.20 

-  1.83 

0.005 

KURTOSIS 

0.61 

0.60 

5.54 

-  0.56 

RANGE 

67.28 

66.35 

125.60 

191 

MINIMUM 

59.52 

59.97 

0.0 

3.75 

.10  QUANTILE 

84.19 

84.25 

72.1 

45.7 

.25 

90.10 

90.13 

85.30 

742 

.50 

96.84 

96.74 

94.26 

102 

.75 

102.96 

102.75 

101.52 

132 

i90 

108.66 

108.40 

107.62 

158 

MAXIMUM 

126.80 

126.32 

125.60 

194 

*   some  CEP  =  0.0 


29 


CASE    III.       P0V7ER   DISTRIBUTION 


4  1 

Y   =    10^    •     (^  -    1) 


MOE  CEP    1  CEP    2  CEP    3  CEP    4 


MEAN 

209.61 

229.33 

100.71 

113 

STD  DEV 

190.15 

257.90 

21.07 

46, 

0 

MSE 

2460 

83000 

444.36 

2279 

GED  MEAN 

177.76 

185.62 

98.60 

103 

VARIANCE 

3615.71 

66515.7 

443.86 

2110 

COEF  VAR 

0.90 

1.12 

0.21 

0. 

40 

SKEWNESS 

6.39 

7.32 

0.84 

0. 

76 

KURTOSIS 

57.54 

72.76 

2.05 

1, 

48 

RANGE 

2264.4 

3220.29 

158.52 

337 

MINIMUM 

74.59 

75.08 

48.32 

17. 

2 

.10  QUANTILE 

108.29 

108.79 

76.88 

57. 

1 

.25 

126.07 

127.96 

86.19 

82. 

5 

.50 

161.63 

166.20 

99.25 

107 

,75 

220.09 

228.03 

112.71 

140 

.90 

331.77 

364.51 

124.77 

172 

MAXIMUM 

2339.01 

3295.38 

206.84 

354 
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CASE  IV.   TRANSLATED  RAYLEIGH  DISTRIBUTION 


Y  =  -  14282.6809  In  U  +  100 


MOE 

CEP  1 

CEP  2 

CEP  3 

CEP  4 

MEAN 

99.79 

99.79 

95.85 

108 

STD  DEV 

11.56 

11.58 

12.64 

38.3 

MSE 

133.74 

134.09 

176.96 

1524 

GED  MEAN 

99.12 

99.11 

94.78 

100 

VARIANCE 

133.70 

134.05 

159.74 

1460 

COEF  VAR 

0.12 

0.12 

0.13 

0.35 

SKEWNESS 

0.15 

0.15 

-  0.56 

0.36 

KURTOSIS 

0.10 

0.08 

4.59 

0.34 

RANGE 

70.75 

70.25 

135.40 

248 

MINIMUM 

67.44 

67.46 

7.73 

20.4 

.10  QUANTILE 

84.81 

84.78 

81.18 

59.0 

.25 

92.04 

92.00 

87.60 

82.7 

.50 

99.48 

99.48 

96.18 

105 

.75 

107.28 

107.12 

103.97 

132 

.90 

114.56 

114.59 

110.99 

156 

MAXIMUM 

138.19 

137.71 

143.13 

268 
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CASE  V.   ^^inULL  PISTRTPUTION   w  -  2 

1 
Y  ■  (-  144260504.  T  In  U)  "^ 


)E 

CEP    1 

CEP    2 

A 

CEP    3 

CEP    4 

:an 

86.18 

85.80 

8^.25 

100 

:d  dev 

5.23 

5.19 

11.47 

19.6 

]?E 

218.34 

228.61 

379.57 

382.00 

':d  mean 

86.03 

85.64 

98.1 

HRIANCE 

27.35 

26.97 

131.51 

382 

<)EF   VAR 

0.06 

0.06 

0.14 

0.20 

5  NEWNESS 

0.029 

0.028 

-      5.47 

-       0.24 

JJRTOSIS 

0.28 

0.27 

38.40 

0.17 

vNGE 

32.67 

32.4  6 

112.79 

121 

!  NIMUM 

71.49 

71.14 

0.0 

420 

10    QUANTIT.E 

79.37 

79.04 

78.33 

74.2 

25 

82.87 

82.52 

81.89 

88.5 

50 

86.29 

85.89 

85.42 

100.0 

75 

89.50 

90.12 

88.69 

113 

90 

92.50 

92.07 

91.56 

124 

r. XI  MUM 

104.16 

103.61 

112.79 

163 
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CASE  VI.   WEIBULL  DISTRIBUTION  w=1.5 


Y  =  (-  1442695.04  In  U) 


MOE  CEP  1  CEP  2 

MEAN  89.58  89.28 

STD  DEV  7.07  7.04 

MSE  158.65  164.46 

GED   MEAN  89.30  89.00 

VARIANCE  50.0  7  4  9.54 

COEF   VAR  0.0  8  0.08 

SKEVVNESS  0.0  7  0.0  7 

KURTOSIS  0.22  0.21 

RANGE  44.55  44.26 

MINIMUM  69.36  69.10 

.lOQUANTILE  80.27  79.99 

.25  85.39  85.12 

.50  89.59  89.25 

.75  94.37  94.07 

.90  98.31  98.00 

MAXIMUM  113.91  113.37 


CEP  3 

CEP 

4 

87.51 

102 

12.85 

25. 

7 

321.29 

664 

98. 

4 

165.09 

660 

0.148 

-  0. 

25 

-  4.27 

-  0. 

04 

27.91 

0. 

11 

125.57 

161 

0.0 

31. 

5 

78.95 

67. 

6 

83.80 

86. 

4 

88.60 

101 

93.44 

120 

97.48 

133 

125.57 

193 

33 


CASE  VII.   WEIBULL  DISTRIBUTION  w=  .5 
y  =  (-  144.27  In  U) ^ 


MOE 

CEP  1 

CEP  2 

CEP  3 

CEP  4 

MEAN 

166.57 

170.05 

87.45 

146 

STD  DEV 

42.61 

44.99 

31.55 

99.2 

MSE 

6246.94 

6930.96 

1153.15 

11956 

GED  MEAN 

161.26 

164.32 

81.72 

113 

VARIANCE 

1815.38 

2023.96 

995.65 

9840 

COEF  VAR 

0.26 

0.26 

0.36 

0.68 

SKEWNESS 

0.62 

0.70 

0.76 

1.39 

KURTOSIS 

0.53 

0.76 

1.48 

3.26 

RANGE 

250.79 

258.53 

218.72 

746 

MINIMUM 

75.81 

78.48 

18.67 

3.24 

.10  QUANT I LE 

116.20 

117.74 

49.30 

41.2 

.25 

135.54 

137.27 

66.56 

75.8 

.50 

163.59 

165.51 

84.41 

126 

.75 

191.58 

196.45 

105.80 

194 

.90 

220.85 

227.60 

127.69 

282 

MAXIMUM 

326.60 

337.01 

237.40 

750 
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CASE  VIII.   WEIBULL  DISTRIBUTION  w=  .1 
Y  =  (-  3.623886  In  U)  -^^ 


MOE 

CEP 

1 

CEP 

2 

CEP  3 

CEP  4 

MEAN 

2.6 

10  = 

3.6 

10^ 

131.06 

5.1  lo"^ 

STD  DEV 

7.8 

10  = 

1.1 

10^ 

335.02 

2.8  10^ 

MSE 

6.7 

10^1 

1.4 

10^2 

1.13  10^ 

8.2  10^0 

GED  MEAN 

6.0 

10^ 

7.4 

10^ 

85.89 

674 

VARIANCE 

6.0 

10^1 

1.3 

10^2 

1.1  10^ 

7.9  10^0 

COEF  VAR 

2.97 

3.14 

2.56 

5.43 

SKEV7NESS 

7.95 

8.38 

14.16 

14.4 

KURTOSIS 

81.28 

89.17 

232.51 

255 

RANGE 

9.8 

10 

1.4 

10^ 

6214.98 

5.4  10^ 

MINIMUM 

461. 

12 

476. 

59 

5.65 

7.2  10^ 

.10  QUANTILE 

685S 

).66 

7661.19 

30.28 

4.41 

.25 

1.8 

10^ 

2.2 

10^ 

54.43 

77.0 

.50 

5.9 

10^ 

7.4 

10^ 

93.37 

1070 

,75 

1.7 

10  = 

2.2 

10^ 

140.28 

9280 

.90 

5.6 

10  = 

7.3 

10^ 

195.25 

8.3  lo"^ 

MAXIMUM 

9.8 

10« 

1.4 

10^ 

6220.63 

5.4  10^ 
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COMPARISON   OF   ESTIMATORS 
USING   MEAN    SQUARED    ERROR  AS    A    CRITERION 

MSE[CEP]    =    S^^p    +     [E(CEP)     -    CEP]^ 


CASE  NO. 

DISTRIBUTION 

CEP  1 

CEP  2 

A- 

CEP  3 

CEP  4 

I 

RAYLEIGH 

136 

136 

165 

1549 

II 

UNIFORM 

109 

107 

368 

1744 

III 

POWER 

2460 

83,000 

444 

2279 

IV 

TRANSLATED 

RAY. 

134 

134 

177 

1524 

V 

WEIBULL 

w 

=  2. 

218 

228 

379 

382 

VI 

WE I BULL 

w 

=  1.5 

158 

164 

321 

664 

VII 

WEIBULL 

w 

-  .5 

6246 

6930 

1153 

11956 

VIII 

WEIBULL 

w 

=  .1 

6.77  lO^-*- 

6.78  1 

0^1 

1.13  1 

.(? 

8.17  10 

10 
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VII.   CONCLUSION 

In  general,  we  found  the  maximum  likelihood  power 
estimator  CEP  3  is  better  than  moment  method  power  estimator, 
in  that  it  does  not  have  the  limitations  and  exists  for  all 
parent  distributions  we  considered.   Using  minimum  mean 
squared  error  as  a  criterion,  the  estimators  based  on  the 
power  model  have  inferior  performance  in  the  "shallow  tailed" 
cases,  but  were  best  in  the  heavy  tailed  cases,  this  is 
especially  true  for  the  w  =  .1  Weibull  case,  in  which  it 
has  less  bias  than  its  competitors  and  has  relatively  small 
variance . 

We  found  the  Rayleigh  based  CEP  estimators  (CEP  1)  to 
be  possibly  the  best  overall,  especially  for  the  shallov; 
tailed  cases.   The  maximum  likelihood  jackknife  Rayleigh 
estimator  (CEP  2)  had  no  better  performance  than  the  unbiased 
Rayleigh  estimator. 

The  value  of  the  sample  median  estimator  (CEP  4)  can 
be  found  for  small  samples  with  such  rapidity  that  in  certain 
cases  the  time  saved  may  compensate  for  the  accuracy  lost. 

In  practice,  when  analysts  find  several  impact  points 
comparatively  far  away  from  target,  they  should  choose 
maximum  likelihood  power  estimator  (CEP  3) .   Otherwise  they 
could  choose  Rayleigh  based  CEP  estimator  (CEP  1) . 
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APPENDIX    l:       MAXIM.UM    LIKELIHOOD    PCWcR    ESTIMiTHR 


10 

3£ 


D  I  ^' ENS  I  ON    Y(  200)  t  S(2001),A(  2001) 

DI^'ef^lSIC■^J  cepcsoo) 
Call  GVFLi^ 

ISEED=470723 

READ    (5,100)     AO,Sn 

RE/D    (5,200)       Ni, IST3P,S(1) ,A(1) 

SN=ELGAT (MS) 

F^=-1.0/A J 

KS=500 

DC     1    KI=1,     KS 

GENER/iTE    Pni^ER    SQUARE    DEVIATE 

DC    5    i\=l,MS 

Y(rvi)  =  0.0 

CALL    RANDOM    (I?EEDtU,l) 

y(,n)  =  s:j--(u-<=^fa-i  .0  » 

CCMTNuE 

SOLVE    NCNLINEAF.     pQUATICN    JNDEk    NEWTON'S     SEC^.^JL    METHOD 

DG    10    I=l,I>TOi' 

ITERS=I 

A(I+1)-0.0 

S(I+1)-C.0 

C^LL    SF     (Si:  ),YtMS,B) 

CALL    CF     (S(  •  ),  Y,.iS,C) 

CALL     DF     (S(I ), Y,NS,0) 

E  =  flLCG( S(  I)  ) 

UX=Si\/A(T  )+3r^^c-i: 

VX  =  Si\;-A  (  I  )/£  (I  )-(  A(I  )  +  1.0)-D 

WX  =  -SN/(A(  T  )-'.'-^^2  ) 

XX  =  Sf!/S  (  I  )-D 

ZX  =  -S-N|-A(  I)/(S{  I  j-^'^?)*!  M  I  )+1.0)--i'C 

DET  =  KX-ZX-XX-^:^2 

CHECK    THE    DETFP MI N&TI ON 
!F     (CET.EQ.;-.0)     GO    TO    ^9 
UU  =  ZX*UX-XX--^\/X 
VV=WX*VX-XX^JX 

WHETHER     TO    STOP    "HE     ITERATION 

CHT=-(  LX-ULl-tVX-VV)/DET 

IF     (CHI .LE.  ?.:4J  )     GO    TO    0  3 

A(  I  +  l  )  =  A(  I  )-'JU/')^T 

S(  I  +  l )=S( I  )-VV/OcT 

IF    (S(I+1 )  .GT.0.0)    GO    TO    10 

S(!-H1)  =  S(  I  )/2.0 

CCNTI'JLIE 

ACU"=A( ITEkS ) 

SGUV=S(  ITE?^S  ) 


ESTIMATE    CEP 
AK  =  2--'^(  1.0 /A  OUT) -1.0 
C  E  P  (  K  I)  =  (  S  OUT  *AK  )  *  -0  .  5 

WRITE     (6,40J)    SOUT, AjUT, ITERS, CEP(K! ) 
1    CONTINUE 

CALL    HISTG    (CEP,K5,0) 
99    STOP 
END 
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SUBROUTINE  FF  (S,Y,N3,3) 

DI^'E^iSION  Y(200) 

B=C.0 

DC  50  N=1,NS 

B=e+ALuG(S+Y(N)) 

CCNTINUE 
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RETURN 
END 

SL6P>CUTI^JE    CF     (StY,N'"^,C) 

DlMEf^SICN    Y(200) 

C  =  0.0 

DO    40    r.^ltNS 

C=C4-1.0/(S+Y(M)  )  ^=*2 

RETURN 
END 

SUBROUTINE    OF    (S,Y,MS,0) 
DIVE.NSICN    Y(200) 
0  =  0.0 

02    50    i\  =  l,NS 
D=D+l.C/(5+Y(N) ) 
50    CTNViMUe 
RETURN 
END 
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APPENDIX    2:    RAYLEIGH    UNPT<iSEO    AND    J£kCKKNIFF    ESTIMATOR 


C 

c 
c 


c 
c 

c 

c 
c 


c 

c 


c 
c 


CI  MENS  I  OM    Y(50,50)  ,SP(  50),  PSE(  50)  »CF.P(  5  00) 
DIME.NSICN    RAY  (300) 

GENERATE    RAYLEIGH    DEVIATE 
CALL    OVFLOW 
ISEED=470723 

f^ead  i\'  original  pgp.  parameter  vc  t  total  savple  siz^ 
nSt  suegrcup  sampli^  size  KS 
read  (5,100)  v.j,ns,ks 

CPC^l.  17741C02--;=SQkT(  VO  ) 
MS=NS/K5 
Sf^  =  FLOAT(MS) 
SN  =  FlO.:..T(NS) 
SK=FLuAT(KS) 

WRITE     (o,2CJ)    VD,CPC,i\'S,MS  ,KS 
KE=5C0 

00    1    KCE=1,KE 
DO    5    M=1,^^.S 
DO    5    K=1,KS 
Y  (  V  ,  K  )  =  0  .  0 

CALL    RANDJVi     (ISEEDtU,  1) 
Y ( V , K ) = 10 000. 0- ( 1 . J/U- 1.0) 
5    CGKTINUE 

ESTIMATE    THE    PARAMETER    BASED    ON    ALL    DATA 
T  =  C.O 

DO    10    N=l,Mr 
CO    10    K=1,KS 
T  =  T+Y(y;,K) 
10    CONTINUE 

TN,  =  SCRT(T/(SN--"=2.0)  ) 

NAXIM'J^-"     LIKELIHOOD    UNBIASED    ESTIMATEP 
RAY(KRE  )=1.184b*TM 

ESTIMATE    TH'r;     PARAMETER    3ASED    ON    ALL    DATA    EXCEPT    THOSE 
IN    THE     ITH    GROUP 

CO   20   y>=i,N*s 
s  =  o.o 

DC    25     1=1, MS 


IF    (I.EQ.M)    GO 


25 


DC    26    ^=1tKS 

S=S+Y( I ,K) 
26    CCKTINJE 
2  5    CONTINUE 

SP(V')^SOkT(S/(  (SN-SK)  '^2.0)  ) 
20    CONTINUE 

COMPUTE  THE  MS  PSEUDG  VALUES 
DO  hC    •^  =  1,'"^. 

P  S  E  {  M. )  =  S  '^  ^  T  -^1  -  (  S  •-!  - 1 .  0  )  -  S  P  (  M  ) 
AO  CONTINUE 

ESTIMATc    JaCKKNIFE    PARAMETER 
ACK=(..0 
CO    4b    y=lTMS 
ACK=ACK+PSE(M) 
45    CONTINUE 

VJACK=ACK/SM 
SCK=0,0 

COMPUTER    THE    CONFIDENCE    INTERVAL 
00    50    1=1, MS 

SCK=SCK+(PSE( I )-VJACK)**2 
50    CCKTINUE 
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SVJA=SCK/(SM=!'(?M-1.0)  ) 
UPE  =  VJ-^CK-2.D^=SQRT  (SVJM 
PLB=VJACK+2.0*SgKT(SVJ4) 

ESTlMATr    C5?    FR?M    J^CKKNI^E     PARAMETER 

CEP(KRr  )=1.1774=^VJACK 

ViRlTE     (6,30JJ     VJACK,UPB,PL6,CEP(KkE) 

CCMIiNUE 

C/^LL    HISTG    (kAY,KE,0) 

C6LL    HISTG    (CEPtKFrO) 

STOP 

END 
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